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ABSTRACT 

Dedicated monitoring of the transient X-ray pulsar XTE J1543-568 during the first 
year after its discovery has revealed the unambiguous detection of a binary orbit. The 
orbital period is 75.56 ± 0.25 d, and the projected semi-major axis 353±8 lt-sec. The 
mass function and position in the pulse period versus orbital period diagram are con- 
sistent with XTE J1543-568 being a Be X-ray binary. The eccentricity of less than 
0.03 (2a) is among the lowest for the 12 Be X-ray binaries whose orbits have now been 
measured. This confirms the suspicion that small kick velocities of neutron stars during 
supernovae are more common than expected. The distance is estimated to be larger 
than 10 kpc, and the luminosity at least 10 37 erg s —1 . 

Subject headings: X-rays: stars — binaries; general — pulsars: individual (XTE J1543- 
568) 

1. Introduction 

The transient X-ray source XTE J1543-568 was discovered during a slew of the Proportional 
Counter Array (PCA) on the Rossi X-ray Timing Explorer (RXTE) on Feb. 6.84, 2000, when 
it had a 2 to 10 keV flux of 2 x 10~ 10 erg cm~ 2 s -1 (Marshall, Takeshima & In 't Zand 2000). 
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A subsequent pointed PCA observation on Feb. 10.69 showed a similar intensity and revealed a 
pulsar with a period of 27.12±0.02 s. Later, the pulsar was found in earlier data from the Burst 
and Transient Source Experiment (BATSE) on the Compton Gamma-Ray Observatory; the first 
detection of the pulsar was on Jan. 25 (Finger & Wilson 2000). The most accurate location was 
recently determined by serendipitous observations with the Wide Field Cameras on BeppoSAX: 
"2000.0 = 15 h 44 m 01 s , £2000.0 = _ 56°45J9 with a 99% confidence region delimited by a circle 
with a radius of 2'0 (Kaptein, In 't Zand & Heise 2001). Thus far, no optical counterpart has been 
reported. 

XTE J1543-568 continues to be active one year after discovery and is being monitored through 
regular pointed PCA observations. We here discuss the results from a timing analysis of data from 
the first year of monitoring, and report the discovery of the binary orbit of XTE J1543-568, present 
the light curve and discuss implications for the companion of the neutron star. 

2. Observations and data analysis 

A target-of-opportunity program was initiated after the discovery with RXTE to monitor the 
pulsar and to try to detect Doppler delays from a binary orbit. The frequency of the observations 
went from once every day for the first 7 days, to once every 2 days during the next 10 days, to 
once every 3 days during the next 75 days, and finally to once every week for the remainder of 
the time. Weekly observations still continued one year after discovery. We here report about the 
first year of observations (i.e., up to Feb. 3, 2001) for which the total accumulated net exposure 
time is 172 ks. The number of observations is 77. Most observations were carried out over two 
consecutive RXTE orbits with a total span of 1.9 hr, but the time span has been as high as 5.5 hr. 
We discuss measurements that were carried out with the PCA on RXTE. This instrument (for 
a detailed description, see Jahoda et al. 1996) consists of an array of 5 co-aligned Proportional 
Counter Units (PCUs) that are sensitive to photons of energy 2 to 60 keV with a total collecting 
area of 6500 cm 2 . The spectral resolution is 18% full-width at half maximum (FWHM) at 6 keV 
and the field of view is 1° FWHM. During the observations the number of active PCUs varied 
between 1 and 5. 



2.1. Light curve 

We extracted spectra for each RXTE orbit of observation, for all active PCUs combined, and 
for layer 1 only. All spectra were successfully modeled between 3 and 20 keV using a power-law 
function with low-energy absorption (following Morrison & McCammon 1983), multiplied with an 
exponential cutoff function and supplemented with a narrow emission line near 6.5 keV which may 
be identified with Fe-K line emission for neutral to C-like ionization states and which has a constant 
equivalent width of 127±5 keV. Figure 1 presents the best-fit spectral parameters and 3-20 keV flux 
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as a function of time, with a time resolution of one observation. The light curve exhibits variability 
on a time scale of weeks. The peak was reached about 90 days after discovery and is nine times 
the lowest flux levels observed thus far. Outside this main peak, the flux occasionally reaches half 
the peak value. There is no obvious periodicity. 

We also studied the light curve as measured with the All-Sky Monitor on RXTE since 1996. 
This light curve is quite noisy because of the small flux. There is no conclusive evidence for 
outbursts other than the one reported here, and there is no indication that XTE J1543-568 was 
ever brighter than the flux levels measured with the PCA. 

2.2. Pulse period and binary orbit 

For the study of the pulse period, we employed event data with 0.5 s time resolution and 
64 readout channels for photon energy (i.e., configuration E_500us_64M_0_ls). We optimized the 
accuracy of the pulse period determination as a function of detector layers and bandpasses and find 
the best result if all layers are included and photons with energies higher than 20 keV excluded, 
although this optimum is only marginally better than if only layer 1 is included. We combined the 
data from all active PCUs, applied standard filters to eliminate bad data, and corrected the resulting 
light curves to the solar system barycenter. The latter includes a correction for the satellite orbit. 
For each of the observations, we measured the pulse period by epoch folding the data for a relevant 
range of test periods and calculating the x 2 values for a constant flux hypothesis. The resulting 
curves show sinc-like trends. To obtain the best-fit period we fitted for each observation the peak 
curve with a 3rd degree polynomial and identified the period with the peak position of the cubic 
function. The 68% confidence error was estimated by determining the range of periods for which 
the x 2 value is within 1 from the maximum. This procedure was verified through Monte Carlo 
simulations. The resulting measurements are presented in figure 2. They are supplemented with 
78 BATSE measurements from the public domain 6 and from Finger & Wilson (2000), pertaining 
to TJD 11572-11678 (TJD = JD-2, 440, 000.5). The BATSE observations stopped prior to de- 
orbiting the satellite in June 2000. Figure 2 clearly shows that there is a 75 d periodicity which 
must be due to an expected binary orbit. Also, long time scale spin-up trends are visible which 
are expected to result from accretion torque. The spin-up is faster in the first ~ 100 days which is 
consistent with an accretion torque explanation because the flux and, therefore, the accretion rate 
is larger at that time. 

In order to disentangle the orbital modulation from the spin up trend, we modeled in an 
initial step the orbital modulation through a simple sinusoidal variation. The spin up was modeled 
through various simple functions with similar success and we here present the results when the 
period derivative P is described by a Gaussian function on top of a flat, non-zero, level. We only 
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modeled the data after TJD 11610 because the spin up seems to be more complex before that. For 
the determination of the orbit this is not detrimental because the excluded time interval represents 
only about 5% of the complete time span of the observations. The results of the modeling are 
presented in figure 2. In a second step, we subtracted the model for the spin up from the data 
and fitted a model of an eccentric binary orbit to the residuals. The best-fit xl is rather high at 
10.2 (for 118 degrees of freedom u), which we attribute to an incomplete spin-up model. In order 
to obtain estimates for the errors in the model parameters that take into account the systematic 
uncertainties, we forced xi to 1 through a multiplication of all pulse period errors with a constant 
factor \/10.2 = 3.2. The results for the orbital and spin- up parameters are given in table 1. In 
figure 2 (bottom panel) we plot the data after subtraction of the orbital modulation model, in 
figure 3 after subtraction of the spin- up model and folding with the orbital period. As a double 
check, we carried out the same procedure for the data after TJD 11670 which excludes the times 
when the Gaussian component of the spin-up model is important. Thus, possible strong systematic 
errors are avoided at the expense of accuracy. The result, given in an extra column in table 1, is 
consistent with the previous result. 

The pulse profile was studied for 3 observations throughout the total observation span, see 
figure 4 for the profile at flux peak. The profile consists of two broad peaks at a phase difference of 
0.5 and with different magnitudes. The pulsed fraction is consistently 65±5%. This kind of profile 
and amplitude is not unusual for pulsars among high-mass X-ray binaries (c.f., Bildsten et al. 1997 
and White et al. 1995). We determined the profiles in two energy bands (2 to 10 keV and 10 to 20 
keV) and do not find strong morphological differences. The pulse fraction is about 10% higher in 
the upper energy band. 



3. Discussion 

Given the projected semi-major axis a x sini and P or b, we can place constraints on the companion 
of the neutron star. The mass function is 

4^ 2 (q x sim) 3 _ (M c sim) 3 m 

If M x = 1.4 M , this implies a minimum mass for the companion star of 10.0 M Q . For a main 
sequence star, the implied spectral type is B3 or earlier (Zombeck 1990). The orbital period places 
XTE J1543-568 in the 'Corbet' diagram of pulse period versus orbital period unambiguously on 
the branch of optically confirmed Be transients (Corbet et al. 1986, c.f. Bildsten et al. 1997). 

In principle the distance is constrained by the spin-up rate and flux. The flux is a measure 
of the mass accretion rate which determines the torque on the neutron star and the spin-up rate. 
Given the steadiness of the spin-up rate, we consider it likely that the accretion occurs through a 
disk. In that case Ghosh & Lamb (1979) predict a spin-up rate of 

-P = l.9x^- l2 ^% 7 m'^Rfl^P 2 L% 7 ss- 1 (2) 
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for a neutron star with magnetic dipole moment ^30 in units of 10 30 G cm 3 , mass m in units of 
1.4 Mq, radius Rq in units of 10 km, moment of inertia I45 in units of 10 45 g cm 2 , and luminosity 
L37 in units of 10 erg s _1 . In figure 5 we present the period derivative as a function of 3-20 keV 
flux. Our data are rather noisy and do not cover a large dynamic range in flux, but they appear 
to be roughly consistent with a 6/7 power law relationship. Given the spectral shapes found, the 
overall bolometric correction is about a factor of 2. If we calibrate the above — P/L37 relation at 
the — P = 1 x 10~ 8 s s~ x level and assume all normalized parameters to be 1, the implied distance 
is 26 kpc. This would place XTE J1543-568 at an unlikely position. Acceptable distances would 
perhaps be at least about 50% smaller. It is hard to achieve that by adjusting the other parameters 
individually: ^30 would have to be at least 64 times larger, m 16 times smaller, Rq 4 times larger, 
or Z45 3 times smaller. Therefore, given the correctness of the Ghosh & Lamb formulation, there is 
at least the suggestion that the distance is large. 

Another constraint for the distance, that does not require fi, may be obtained from an argument 
over the magnetospheric radius r mag at which the neutron star magnetosphere begins to dictate the 
mass flow. This radius scales with M -2 / 7 or F 2 / 7 . If the flux F varies between F m j n and F max at 
times when the mass accretion persists, r mag < r co (-F m i n /.F max ) 2 / 7 at peak flux, where r co is the co- 
rotation radius at which the Keplerian period is equal to the spin period (r co = (GMP 2 /47r 2 ) 1 / 3 ). 
The torque 2n 1 v is given by M ^GMr ma _ g . Substition in the above condition reveals 

M > -5.6 x 10 2 ——-I^- i _ M yr" 1 (3) 

m 2/3p7/3 (Fm . n/jPmax) l/7 O k J 

For standard parameters and F max = 4.5F min (after TJD 11670), we obtain M > 9.8x 10~ 10 M yr _1 . 
If all liberated potential energy is released in the form of radiation, the bolometric luminosity is 
constrained by L > 1.2 x 10 37 erg s _1 . For F m3iX = 4.5 x 10~ 10 erg s _1 cm~ 2 and a bolometric 
correction of 2, a minimum distance of 10 kpc is implied. Again this is a fairly large distance. 

These distance estimates suggest that searches for the optical counterpart may prove difficult, 
also because the large value of A^h points to a high extinction. For a 10.0 M main sequence star, 
My = —2. If Nr = 1-3 x 10 22 (see figure 1; this is consistent with the total Galactic column as 
inferred from maps by Dickey & Lockman 1990), Ay = 7.4 (according to the standard conversions 
introduced by Predehl & Schmitt 1995). For a distance of 26 kpc, the apparent visual magnitude 
would be V = 23, for a distance of 10 kpc, V = 21. 

The most interesting result of our analysis is the low eccentricity. It is among the lowest for 
the 12 confirmed or suspected Be X-ray binaries whose orbits have now been determined (8 cases 
are reviewed by Bildsten et al. 1997 and since then orbits have been determined for 2S 1845-024 
by Finger et al. 1999, for SAX J2103. 5+4545 by Baykal, Stark & Swank 2000, and for X Per by 
Delgado-Marti et al. 2001). There is one other system apart from XTE J1543-568 which has an 
eccentricity that is consistent with (e < 0.09 at 2a confidence for 2S 1553-54, see Kelley et al. 
1983). However, one should note that Kelley et al. have a sparse coverage of less than one binary 
orbit and do not solve for spin-up. Finding low-eccentricity systems is remarkable because one 
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would expect most long-period Be systems to have undergone supernovae with large kick velocities 
for the neutron stars and to have circularization time scales far longer than their age. XTE J1543- 
568 confirms the suspicion that there may be a substantial population of neutron stars formed with 
little or no kick (Delgado-Martf et al. 2001). 

We thank the project scientist Jean Swank for continuing the RXTE observations of XTE J1543- 
568 through public TOOs after our TOO program time was spent, Evan Smith for his tireless efforts 
to schedule the many observations, Mike Tripicco for his assistance in the barycentric corrections for 
2001 PCA data, Ron Remillard for his advice on the ASM light curve, and the referee Mark Finger 
for valuable suggestions. JZ acknowledges financial support from the Netherlands Organization for 
Scientific Research (NWO). 
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Table 1: Fit parameters to pulse period history. The errors are for 68% confidence. The numbers 
in the third column reference, from data after TJD 11670 (40 data points) 



Orbital period _P or b (days) 
Orbital amplitude (ms) 
Semi-major axis a x sin* (lt-sec) 
Eccentricity e 

Epoch of mean longitude tt/2 (TJD) 
Pulse period P at TJD 11600 (s) 
Constant P term (s s _1 ) 
Gaussian P term 
a (days) 

peak time (TJD) 
peak P value (s s _1 ) 



75.56±0.25 

9.18±0.21 

353±8 

< 0.03 at 2a 

11725.88±0.16 

27.12156±0.00059 

-(2.84 ±0.02) x 10- 9 

13.30±0.02 
11636. 71±0.02 
-(1.27 ±0.03) x 10~ 8 



75.88±0.40 

9.07±0.26 

349±20 

< 0.03 at 2a 

11725. 54±0.62 

27.1225±0.0043 

-(2.84 ±0.03) x 10~ 9 
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Fig. 1— Top panel: PCA-measured time history of XTE J1543-568 flux in 3 to 20 keV, from 
the time of discovery through the first year of monitoring observation. The errors (not indicated) 
are typically 1%. Second through sixth panel: evolution of power-law photon index, e-folding 
cutoff energy, photon flux and centroid energy of the emission line, and hydrogen column density 
For some spectra, the cutoff energy could not be constrained and no data point is given in the 
appropriate panel. The average hydrogen column density is (1.3 ± 0.1) x 10 22 cm 2 . The emission 
line flux follows the overall flux. In fact, the equivalent width is consistent with being constant 
with a value of 127 eV. Bottom panel: % 2 statistic for the relevant spectral fits. 
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Fig. 2. — Time history of pulse period as measured with the PCA and BATSE. The solid curve 
shows the fitted model (see text) that was applied to data after TJD 11610. The middle panel 
shows the residuals of the data points with respect to the fitted model in terms of a. The bottom 
panel shows the pulse period history after subtraction of the sinusoidal variation. All error bars 
are without multiplication by 3.2 (see text). 
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Fig. 3. — Pulse period history for data taken a 11610, after subtraction of spin-up model 

and folding with period of 75.56 days. The solid curve shows the solution for the binary orbit, with 
e = 0.005. All error bars are without multiplication by 3.2 (see text). 
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Fig. 4. — Pulse profiles for the observation taken on TJD 11637, when the peak flux was reached. 
The solid histogram is for 3-10 keV, the dashed one for 10-20 keV. 
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3 = 20 keV flux [erg s"' cm" 2 ] 

Fig. 5. — Period derivative, as determined from BATSE and RXTE data points, versus 3-20 keV 
flux. The baseline for each period derivative determination was three observations (or at least 2 
days). The error bars on the flux represent the range of fluxes measured during this baseline. The 
solid line is the best fit power-law function, with index 0.95 ± 0.06. 



